Transition Metal-Ethylene Complexes as High-Capacity Hydrogen Storage Media 
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From first-principles calculations, we predict that a single ethylene molecule can form a stable 
complex with two transition metals (TM) such as Ti. The resulting TM-ethylene complex then 
absorbs up to ten hydrogen molecules, reaching to gravimetric storage capacity of ~14 wt%. Dimer- 
ization, polymerizations and incorporation of the TM-ethylene complexes in nanoporous carbon 
materials have been also discussed. Our results are quite remarkable and open a new approach to 
high-capacity hydrogen storage materials discovery. 

PACS numbers: 73.63.Nm, 72.25.-b, 75.75.-|-a 
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Hydrogen is considered as one of the best alternative 
and renewable fuelsfH, ^ due its abundance, easy syn- 
thesis, and non-polluting nature when used in fuel cells. 
However, the main concern is the safe storage and effi- 
cient transport of this highly flammable gasj^l- 

The main obstacles in hydrogen storage are slow kinet- 
ics, poor reversibility and high dehydrogenation temper- 
atures for the chemical hydrides'^; and very low desorp- 
tion temperatures/energies for the physisorption materi- 
als (metal-organic frameworks (M0F)[3, carbide-derived 
carbons (CDC)@|, etc). Recently, a novel concept to 
overcome these obstacles have been sugg ested[l|i|i[ll 
llll ll^ . It was predicted that a single Ti-atom 

affixed to carbon nanostructures, like Cgo or nanotubes, 
strongly adsorbs up to four hydrogen molecules (7t. i8t. ilQi] . 
The interaction between hydrogen molecules and tran- 
sition metals is very unique, lying between chemi and 
physisorption, with a binding energy of 0.4 eV compati- 
ble with room temperature desorption/absorption. The 
origin of this unusual "molecular chemisorption" is ex- 
plained by well-known Dcwar coordination and Kubas 
interaction 15]. The transition metals (TM) are chem- 
ically bonded onto different moleculcs/nanostructures 
through hybridization of lowest-unoccupied molecular or- 
bital (LUMO) of nanostructure with TM c?-orbitals (i.e. 
Dewar coordination). The resulting complex then binds 
multiple molecular hydrogens through hybridization be- 
tween H2 — (7* antibonding and TM d-orbitals (i.e. Kubas 
interactions). 

Synthesizing the predicted structures of Ti decorated 
nanotubes/Ceo were proven to be very difficult due to 
lack of bulk quantities of small-diameter nanotubes and 
strong Ceo-Cgo interactions in the solid phase. Moreover, 
Ti atoms uniformly coating the SWNT/Cgo surface may 
be subject to clustering after several charging-discharging 
process 0. 

In search for a more efficient and feasible high ca- 
pacity hydrogen storage medium, we found that the 
C=C double bond of an ethylene molecule C2H4, mim- 




FIG. 1: (color online) (a) One of the most stable structure 
of Ti-Ceo complex where Ti atom (blue) is bonded to a dou- 
ble bond with four hydrogen molecules attached (red), (b) 
The local structure of the Ti-Ceo-double bond, (c) Replacing 
the end carbon atoms shown in (b) by H results in an ethy- 
lene molecule. This suggests that we may simple use ethylene 
molecule to hold Ti atoms, which then binds multiple hydro- 
gen molecules. 



ics double bonds of Ceo that strongly binds TM-atom 
(see Fig. 1) and therefore it is expected to support TM 
atoms strongly to provide a basis for high-capacity hydro- 
gen storage via the Dewar-Kubas mechanism discussed 
above. 

In this letter, we explored this idea and indeed found 
that a single ethylene molecule can hold not only one 
but two Ti atoms, i.e. C2H4Ti2, which then reversibly 
binds up to ten H2 molecules yielding an unexpectedly 
high storage capacity of 14 wt %. These results suggest 
that ethylene, a well-known inexpensive molecule, can be 
an important basis in developing frameworks for efficient 
and safe hydrogen storage media. 

Our results are obtained from first-principles plane 
wave calculations within density functional theory using 
Vanderbilt-type ultra-soft pseudopotentials with Perdew- 
Zunger exchange correlation ^17^. Single molecules have 
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FIG. 2: (Color online) Optimized structure of an ethylene 
molecule C2H4 (a), C2H4Ti (b), and C2H4Ti2 (c). The panel 
(d) shows that Ti-C2H4 bonding orbital results from the hy- 
bridization of the LUMO of the C2H4 and the Ti-d orbital, in 
accord with Dewar coordination. 

been treated in a supercell ofl5xl5xl5A with T k-point 
and a cutoff energy of 408 eV. The structures are opti- 
mized until the maximum force allowed on each atom is 
less than 0.01 eV/A for both spin-paired and spin-relaxed 
cases. 

We first studied the bonding of a single Ti atom to an 
ethylene molecule to form C2H4Ti (see Fig. 2(b)). We 
found no energy barrier for this reaction. The binding 
energy is calculated by subtracting the equilibrium to- 
tal energy Et of C2H4Ti molecule from the sum of the 
total energies of free ethylene molecule and of Ti atom; 
EB(Ti)=ET(C2H4)-hET(Ti)-ET(C2H4Ti). The Ti atom 
forms a symmetric bridge bond with the C=C bond of 
ethylene with Eb =1.45 eV. Interestingly, it is also pos- 
sible to attach a second Ti atom to the C2H4Ti to form 
C2H4Ti2 (see Fig. 2(c)) without any potential barrier 
and about the same binding energy as the first Ti atom. 
In the optimized structure (Fig. 2(c)), each Ti atom is 
closer to one of the carbon atoms, leading to two differ- 
ent Ti-C bonds. Fig. 2(d) shows that the bonding orbital 
for the Ti-atoms and C2H4 results from the hybridiza- 
tion of the LUMO orbital of the ethylene molecule and 
the Ti-d orbitals, in accord with Dewar coordination^^. 
The spin-polarized calculation gives 1.53 eV lower energy 
than non-spin-polarized one, suggesting magnetic ground 
state for C2H4Ti2 with moment fi — 6fiB per molecule. 

The stability of the C2H4Ti2 complex were further 
tested by normal mode analysis. We found no soft (i.e 
negative) mode. There are three main Ti-modes. In two 
of these modes, Ti atoms vibrate parallel and perpen- 
dicular to the C=C bond and have the energies of 176 
cm~^ and 123 cm~^, respectively. In the third Ti-mode 



Ti atoms vibrate perpendicular to the C2H4 plane with 
energy of 367 cm~^. These three modes are unique for 
the C2H4Ti2 complex and therefore should be present 
in any Raman/IR spectra of a successfully synthesized 
material. 

We next studied the H2 storage capacity of Ti-ethylene 
complex, by calculating the interaction between C2H4Ti2 
and different number of H2 molecules and configurations. 
The first H2 molecule is absorbed dissociatively to form 
C2H4(TiIl2)2 as shown in Fig. 3(a) with a binding energy 
of 1.18 eV/H2. The additional hydrogen molecules do 
not dissociate and molecularly absorbed around the Ti- 
atom. Two of these configurations are shown in Fig. 3(b) 
and (c). In C2H4(TiIl2-2Il2)2 configuration, two H2 are 
molecularly bonded from left and right side of the Till2 
group with a binding energy of 0.38 eV/H2 and signif- 
icantly elongated H-H bond length of 0.81 A. It is also 
energetically favorable to add a third H2 molecule from 
the top of the TiH2 group, with a binding energy of 0.4 
eV and bond length of 0.82 A. The resulting structure, 
C2H4(TiIl2-3H2)2, is shown in Fig. 3(c). We note that 
these binding energies have the right order of magni- 
tude for room temperature storage. Since the hydro- 
gens are absorbed molecularly, we also expect fast ab- 
sorption/desorption kinetics. 

Finally, we also observed many other local stable 
configurations where all of the hydrogen molecules are 
bonded molecularly. One of such configuration, denoted 
as C2H4(Ti-5H2)2, is shown in Fig. 3(d). Here the H2 
molecules stay intact and benefit equally from bonding 
with the Ti atom. The total ten hydrogen molecules ab- 
sorbed by a single Ti-ethylene complex, C2H4(Ti-5H2)2, 
corresponds to a ~ 14 wt % gravimetric density, which 
is more than twice of the criterion set for efficient hydro- 
gen storage. By removing the top H2 molecule, we find 
that C2H4(Ti-4H2)2 configuration is also a local mini- 
mum and have a slightly (0.08 eV/H2) higher energy than 
the C2H4(TiIl2-3 112)2 configuration shown in Fig. 3(c). 
Our MD simulations indicate that the system oscillates 
between these two configurations. This is consistent with 
the observation that in Kubas compounds |15|. the dihy- 
drogen (i.e 211) and molecular (i e H2 ) bonding are usu- 
ally found to be in resonance|15j. 

The top hydrogen molecule in C2H4(Ti-5Il2)2, has the 
weakest bonding in the system with Eb=0.29 eV while 
the side H2 molecules have the strongest bonding with 
Eb=0.49 eV/Il2 and significantly elongated H-H bond 
distance of 0.85 A. This suggests the presence of two dif- 
ferent H2-C2H4Ti2 bonding orbitals as shown in Fig 3(e) 
and (f). The first one is the hybridization between the 
top-H2 cr*-antibonding and the Ti-d orbital. The sec- 
ond one is the simultaneous hybridization of the side H2 
cr*-antibonding orbitals with the Ti-d orbital. Since the 
bonding orbitals are mainly between metal d— and hy- 
drogen (T*-antibonding orbitals, the mechanism of this 
interesting interaction can be explained by the Kubas 
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FIG. 3: (Color online) Atomic configurations of an ethylene 
molecule functionalized by two Ti atoms, holding (a) two H2 
molecules which are dissociated, (b) six H2 molecules, and (c) 
eight H2 molecules. Panel (d) shows a configuration where 
ten H2 are bonded all molecularly. The spin-polarized calcu- 
lations gave lower energies (in eV) by 1.5, 0.37,0.16, and 0.06 
for configurations shown in (a)-(d), respectively, suggesting 
magnetic ground state for all cases with moment of « 2fiB- 
The panel (e) and (f) shows the bonding orbital for the top 
(e) and side hydrogen molecules, respectively. Note that the 
hydrogen cr*-antibonding orbitals are hybridized with Ti-d- 
orbitals, suggesting Kubas interaction for the H2-Ti bonding. 



iiiteraction[l5j. 

We also calculated the normal modes of C2H4(Ti- 
5132)2 and did not find any soft modes, indicating that 
the system indeed corresponds to a local minimum. 
Among many vibrational modes, we note that the H2 
stretching mode is around 2700-3000 cm~^ for side H2 
and around 3300 cm~^ for top H2 molecules, significantly 
lower than the 4400 cm^^ for free H2 molecule. Such a 
shift in the mode frequency would be the key feature that 
can be probed by Raman/IR measurement to confirm a 
successful synthesis of the structures predicted here. 

Finally, the stability of C2H4(Ti-5H2)2 structure has 
been further studied by carrying out ab-initio molecu- 
lar dynamic calculations performed at 300 K and 800 K. 
While all ten H2 molecules have remained bound to 
C2H4Ti2 molecule at 300 K, they start to desorb above 
300 K and all of them are desorbed already at 800 K, 
leaving behind a stable C2H4Ti2 molecule. This predic- 
tion suggests that all the stored hydrogen molecules can 
be discharged easily through heating. 

We next discuss the possibility of the dimerization and 
polymerization in the course of recycling and its effect 
on the hydrogen storage capacity. We found that two 
molecules can form a dimer through a Ti-Ti bond as 
shown in Fig. \^sl). The dimer formation is exother- 
mic with an energy gain of E£)=2.28 eV, and leads to a 
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FIG. 4: (color online) (a) Atomic configuration and charge 
density plots of a dimer of C2H4Ti2 molecules linked by two 
Ti atoms one from either molecule. Large, medium and small 
balls indicate titanium, carbon and hydrogen atoms, respec- 
tively. Charge density plots (from left to right) correspond 
to HOMO spin-up, and two spin-down states, (b) Atomic 
configuration of C2ll4Ti2 dimer holding 14 H2 molecules, (c) 
Polymer of C2H4Ti2 molecule with polymerization energy 
Ep. (d) Atomic configuration of the polymer holding 6 H2 
molecule, (e) A C2H4Ti2 molecule adsorbed above the center 
of a hexagon in a (3x3) cell holding 7 H2. (f) Two Ti-ethylene 
complexes adsorbed above and below the center of the same 
hexagon holding 14 H2. 

stable structure. The ground state is ferromagnetic with 
M =6/is. While each Ti atoms at both ends of the dimer 
can bind 5 H2 molecules, two linking Ti atoms can ab- 
sorb only 4 H2 totaling to 14 H2 per dimer. As a result, 
the gravimetric density obtained by the dimer is lowered 
to 10 wt%. The total energy can be further lowered 
by adding more C2H4Ti2 molecules to the dimer, and 
eventually by forming a paramagnetic polymer as shown 
in Fig. EJc). The polymerization energy per molecule 
is calculated to be Ep— 2.69 eV and the H2 storage ca- 
pacity is further lowered to 6.1 wt%. Polymerization did 
not change the binding energy of the H2 on Ti-atom sig- 
nificantly and therefore should not affect the desorption 
temperature. 

In order to prevent C2H4Ti2 molecules to form a pos- 
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TABLE I: The binding energies (in eV) with respect to atomic and bulk energies of various metals (M). The last two rows 
indicate the maximum number of H[2 molecule bonded to each metal and its average binding energy (in eV). 



Property/M 


Sc 


Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


Zr 


Mo 


W 


Pd 


Pt 


Eb (M-atomic) 


1.39 


1.47 


1.27 


0.05 


0.37 


0.83 


1.30 


0.70 


1.41 


none 


1.69 


0.37 


1.18 


1.56 


1.78 


Es(M-bulk) 


-2.72 


-3.66 


-4.13 


-3.57 


-3.20 


-1.74 


-2.53 


-2.19 


-2.25 




-4.44 


-5.84 


-7.18 


-2.24 


-3.56 


max H2/M 


5 


5 


5 


5 


5 


5 


3 


2 


2 




5 


5 


5 


2 


2 


Es(per H2) 


0.39 


0.45 


0.43 


0.35 


0.34 


0.26 


0.41 


0.87 


0.14 




0.57 


0.77 


0.90 


0.58 


0.95 



sible polymer phase during recycling, one can imag- 
ine to incorporate these TM-ethylene complexes in a 
nanoporous materials such as MOf|^ and carbide- 
derived carbons 6]. As an example, we consider a 
C2H4Ti2 adsorbed above and below the center of a 
hexagon of a graphene layer. Here the graphene is taken 
as a prototype system which represents the internal struc- 
ture of a carbon-based nanoporous material. We find 
that TM-ethylene complex can form stable structure with 
graphene surface. As shown in Fig. 0] (e) and (f), sin- 
gle and double €2114X12 molecules assembled on a 3x3 
graphene layer can hold 7 and 14 H2 with average binding 
energy of 0.43 and 0.41 eV/H2, respectively. The bind- 
ing energy of the C2H4Ti2 molecule is found to be ~2 
eV. The actual binding energy in nanoporous materials 
could be even higher due to curvature effects 18]. The 
maximum gravimetric density achieved in this present 
framework is 6.1 wt% 

It is important to know if the results reported above 
for M=Ti hold for other metals. Therefore we are cur- 
rently studying a large number of metals and the details 
will be published elsewhere [l^. Our initial results are 
summarized in Table I, which clearly indicates that most 
of the light TM atoms can be bound to ethylene and each 
of them can absorb up to 5 H2 molecules. Scandium is 
the ideal case but for practical reasons Ti is the the best 
choice of elements. Cr binds very weakly while Zn does 
not bind at all to the C2H4 molecule. Interestingly, Zr 
forms stronger bonding with C2II4 than Ti and can ab- 
sorb up to ten II2 molecularly with an average binding 
energy of 0.6 eV. Heavier metals such as Pd and Pt can 
also form complexes with C2H4 but bind less hydrogen 
molecules with significantly stronger binding energy than 
Ti. Table I also gives the binding energies with respect to 
bulk metal energies. The negative value for indicates 
endothermic reaction. Due to very low vapor pressure 
of many metals, it is probably better to use some metal- 
precursor to synthesize the structures predicted here. 

In conclusion, we showed that an individual ethylene 
molecule functionalized by two light transition metals 
can bind up to ten hydrogen molecules via Dewar-Kubas 
interaction, reaching a gravimetric density as high as 
~14 wt%. We propose to incorporate TM-ethylene com- 
plex into carbon-based nanoporous materials to avoid 



dimerization/polymerization during recycling. Our re- 
sults open a new approach to high-capacity hydrogen 
storage materials discovery by functionalizing small or- 
ganic molecules with light transition metals. 
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